NUMERICAL STUDY OF HALO CONCENTRATIONS IN DARK-ENERGY COSMOLOGIES 
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ABSTRACT 

We study the concentration parameters, their mass dependence and redshift evolution, of dark-matter halos in 
different dark-energy cosmologies with constant and time-variable equation of state, and compare them with "stan- 
dard" ACDM and OCDM models. We find that previously proposed algorithms for predicting halo concentrations 
can be well adapted to dark-energy models. When centred on the analytically expected values, halo concentrations 
show a log-normal distribution with a uniform standard deviation of ~ 0.2. The dependence of averaged halo con- 
centrations on mass and redshift permits a simple fit of the form (1 +z)c = co (M/Mq)", with a ~ —0.1 throughout. 
We find that the cluster concentration depends on the dark energy equation of state at the cluster formation redshift 
Zcoii through the linear growth factor (zcoii ) . As a simple correction accounting for dark-energy cosmologies, we 
propose scaUng cq from ACDM with the ratio of Unear growth factors, co — > cq D+ (zcoii ) / i'+.ACDM (zcoii ) ■ 

a^ 
o 

The properties of dark-matter halos have been shown in many 
— , earlier studies to depend in characteristic ways on their forma- 
i-C tion history. The dark-matter density profile consistently found 
O-dn numerical simulations by (Navarro et al., 1996, 1997) has a 
Q radial scale, rs, which is smaller than the virial radius ryir by the 
5^ concentration parameter c — r^ir/rs- Although there is an ongo- 

ing discussion as to the exact density-profile slope in halo cores, 
d there is agreement that the density profile steepens substantially 
^* going outward across the scale radius. 
. !^ The central density in halo cores was found to reflect the mean 

density of the Universe at the time of halo formation. Since 
^ halos of increasing mass form at increasingly late cosmic epochs 

in hierarchical models of structure formation, the concentration 

parameter of halos at fixed redshift decreases with mass. For 

halos of fixed mass, the concentration parameter increases with 

decreasing redshift because the background density drops. 
Structure growth in the Universe is thus reflected by the core 

densities, or concentrations, of dark-matter halos. Haloes thus 

establish a connection between principally measurable quanti- 
ties, like their central densities, and the cosmological framework 

model because the latter dictates how structures form. 

We have shown in an earlier paper (Bartelmann et al., 2002) 

that halos of given mass and redshift are predicted to have higher 

concentrations in dark-energy cosmologies than in models with 

cosmological constant because they tend to form earlier and thus 

acquire a higher core density. For that study, we used analytic 

algorithms for predicting halo properties which were assumed, 

but never shown, to hold true in dark-energy cosmologies. In 

this paper, we describe numerical simulations undertaken with 

the specific aim of studying how halo properties change in dark- 
energy models, and whether the analytic algorithms proposed 

for other cosmologies can be adopted or adapted to dark-energy 

models. 

Klypin et al. iKlvpin et all 120031) also studied the properties 



'present address: ITA, Universitat Heidelberg, Tiergartenstr. 15, D-69121 
Heidelberg, Germany 



of clusters simulated in the framework of cosmological models 
with dynamical dark energy. However, their cosmological mod- 
els differ from ours for the choice of many parameters. More- 
over, we will analyze here a larger sample of individually sim- 
ulated clusters, covering a completely different range of halo 
masses. 

We describe the dark-energy models used in Sect. 2 and our 
numerical simulations in Sect. 3. Sect. 4 details the determina- 
tion of halo concentrations, their comparison with analytic ex- 
pectations, and their statistical properties. We identify an inter- 
esting dependence between average halo concentrations and the 
density of dark energy at halo formation in Sect. 5, and conclude 
with a summary in Sect. 6. 

2. Dark-energy models 

In this paper, we study and compare results obtained for the 
following cosmological models: an open Cold Dark Matter 
(OCDM) and four flat dark-energy cosmogonies. The latter are 
a cosmological constant (ACDM) model, a dark-energy model 
with constant equation of state (DECDM), and two quintessence 
models, one with inverse po wer-law Ratra-Peebles potential 
(RP, see lPeebles & Ratr'aLl200 2. and references therein) and one 
with SUGRA potential (SUGRA, see Brax & Martin, 2000, and 
references therein). 

In all cases, the matter density parameter today is Q.o = 0.3. 
In the flat cosmologies, the remaining 70% of the critical den- 
sity is assigned to the dark energy at present. The remaining 
cosmological parameters are h — 0.7, Q.i,h^ = 0.022, a Gaussian 
density fluctuations with scale-invariant power spectrum, and no 
gravitational waves. 

Particularly important for us is the normalisation of the per- 
turbation power spectrum, which we set by defining the rms 
density fluctuation level within spheres of 8/i^^Mpc radius, Og. 
In this respect, we follow two approaches, which are normal- 
ising the perturbation amplitude either on large scales with the 
obse rved Cosmic Microw ave Background (CMB) anisotropics 
('e.g. lBennettetal.Ll2(K)l . or on small scales using the observed 



cluster abundance. In the second case, we choose Og = 0.9 in 
all the models, while the ag derived from the CMB is generally 
slightly smaller because of the Integrated Sachs Wolfe (ISW) 
effect affecting the large scale CMB anisotropics in the cos- 
mologies we consider (see Bartelmann et al., 2002, and refer- 
ences therein): we take 0.86, 0.82 and 0.76 for DECDM, RP and 
SUGRA, respectively. These numbers as well as all the inputs 
from the linear evolution of cosmological perturbations to the 
A^-body procedure described later are computed using our dark 
energy oriented cosmologi cal code (Perrotta & Baccigalupi, 
11999), based on CMBfast ISehak & Zaldarr iaga. 1996). ' 

We briefly describe now the dark energy cosmologies we will 
adopt. The equation of state w is a key parameter, describing 
the ratio between the dark energy pressure and energy density, 
and it must be negative in order to effectuate cosmic acceler- 
ation today. In the limit of w being constant, the dark energy 
density evolves in redshift proportionally to (1 That 
is the case for both the ACDM model, whose w = — 1 yields a 
constant dark energy density, and the DECDM model, which we 
choose to have w = —0.6. Note that the latter case is an effective 
scenario in the sense that it is obtained simply by replacing the 
cosmological-constant equation of state, assuming it to be con- 
stant and neglecting any other dark energy features such as its 
cosmological perturbations. Also, it is useful to note that from 
the point of view of the Friedmann equation, the curvature term 
in the OCDM model behaves as a dark energy component with 
constant w — —1/3 (see Bartelmann et al., 2002, and references 
therein). 

The dark energy is consistently described by means of the 
quintessence scalar field (|). The dynamics of its unperturbed 
value and of the linear fluctuation 5(|) obey the Klein-Gordon 
equation 



□(j) + y' = o , 5(n(j)) + y"5(j) = o^ 



(1) 



where □ is the d' Alambert operator in a Friedmann-Robertson- 
Walker (FRW) cosmology, and V{(^) is the quintessence poten- 
tial. 

The RP and SUGRA potentials are given by 



Vrp = 



VsUGRA = 



respectively. The exponential in Vsugra comes from super- 
gravity corrections ( Brax & Martin, 2000) and induces large 
changes in the trajectories with respect to the RP case, as we 
explain now. 

The attractor properties of the tracking trajectories in the 
quintessence scenarios allow to reach the present field value (|)o, 
of the order of the Planck mass Mpianck. starting from a wide set 
of initial conditions for (|) and (j), with the only relevant condition 
that (|),- < Mpianck- For both RP and SUGRA, the ti-acking ti-a- 
jectories are well defined until the quintessence is subdominant 
compared to the matter density, yielding a constant equation of 
state obeying the simple relation 



(3) 



2 + aQ 



In the RP scenario, if we require the pres ent equation of state 
WQ to respect the current constraints (see i Bennett et all I2003L 
and references therein), say wq < —0.8, the exponent must be 
in the range < ag < 0.5, yielding a shallow potential shape. 
At present, the tracking regime is abandoned, but the shallow 
potential does not permit wq to be far from the tracking solution 
of Eq. (|3jl, the deviation being at the 10% level. 
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Fig. 1 . — Equation of state for the RP and SUGRA quintessence models 
studied in this work. 



The SUGRA exponential correction flattens the potential 
shape noticeably at (|) ~ Mpianck. i-C- at the end of the tracking 
trajectory. This means that a given equation of state at present is 
obtained for noticeably higher values of ag than for RP, i.e. the 
dark energy dynamics and thus the cosmological expansion rate 
as a function of redshift are generally much different in the two 
cases. 

Thus, we consider two quintessence models here having the 
same equation of state at the present epoch, but markedly differ- 
ent ones in the past. Specifically, they are tracking quintessence 
scenarios obtained out of the RP and SUGRA potentials, both 
having wq — —0.83 today, with ag = —0.6 and ag = —6.7, re- 
spectively. According to Eq. (|3}, this yields wsugra — —0.23 
and wrp ~ —0.77 in the tracking regime, well representing the 
diversity that dark energy cosmologies can have in the past, even 
if they reproduce the same amount of cosmic acceleration today. 
The redshift behaviour of w in the two cases is illustrated in 
Fig.O 

Note that if w depends on z, the dark energy density evolves 



exp (47lG(j)^) , (2) as 



Pde(z) Pde.o exp 



3 / dz 
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which is easily verified solving the energy conservation equation 
p-f 3//(p+p) = 0. 

The cosmological scenarios we study here are mutually very 
different. The differences are relevant even between the flat dark 
energy models: for DECDM, the equation of state is far from 
— 1, while RP and SUGRA have the same equation of state to- 
day, but differ strongly in the past, as we emphasised above. 
These pronounced differences allow us to recognise their im- 
print in the properties of the dark non-linear halos which we 
study next. In this respect, we obtain the most interesting results 
when the cosmologies above are compared by keeping every- 
thing else fixed, including Og. 



Numerical simulations 



3.1. The simulation code 



We adapted the cosmological code GADGET (Springel et alj 
2001c) for simulating a set of 17 galaxy clusters in various dark- 
energy cosmologies with constant and time-variable equations 
of state. This section describes the code used, the modifications 



required for dark-energy models, and the initial conditions we 
set up. 

In total we performed 136 cluster simulations. The simula- 
tions were performed in parallel on 8 CPUs on an IBM-SP4 
located at the CINECA Super Computing Centre in Bologna, 
Italy. 

3.1.1. GADGET 

The cosmological code GADGET (Springel et al. 2001b) is well 
tested and used for simulating a wide ran ge of cos r nolo gical 
scenarios from the interaction of galaxies ('Soringel', '2000) to 
large-scale structure formation (Springel et al., 2001a). GAD- 
GET is based on a tree-SPH code using co-moving coordinates. 
The new version (P-GADGET-2, kindly provided by Volker 
Springel) allows the computation of long-range forces with a 
particle-mesh (PM) algorithm, with the tree algorithm supply- 
ing short-range gravitational interactions only. This "TreePM" 
method can substantially speed up the computation while main- 
taining the large dynamic range and flexibility of the tree algo- 
rithm. It also differs in the choice of internal variables for time 
and velocity, in its time-stepping algorithm and in its parallelisa- 
tion strategy (see Springel 2003, in preparation). The modifica- 
tions necessary for performing simulations within cosmologies 
with dark-energy within P-GADGET-2 do not differ in principle 
from the modifications necessary for the earlier version of GAD- 
GET, but appear in slightly different places due to the changes 
in design of the code. 

3.1.2. Including dark-energy in GADGET 

Within the co-moving coordinate scheme of GADGET, the only 
place where the dark energy has to be taken care of is the calcu- 
lation of the Hubble function, or some combination of the Hub- 
ble function with the scale factor a. This is used when ever a 
conversion to physical quantities is needed, like converting the 
internal time variable log a to physical time t, or in the equation 
of motion in a cosmological context. Thus, for running cosmo- 
logical simulations including dark-energy, once can rewrite the 
usual Hubble function 



H{a)=Ho 
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1 — — 



1/2 



(5) 



for a flat cosmology with cosmological constant using Eq. (0} as 

1 — ilo — 



Hia) = Ho 



^3 



l+w{a') ^ , 
Q.Q exp I -3 / ; da 



(6) 



Within the new integration scheme of P-GADGET-2, this is also 
done within the calculation of the drift and kick factors (Springel 
et al., in preparation). 

For reasons of timing and precision, we tabulate the integral 
in Eq. (|6} at the beginning of a run, and then interpolate within 
the table during the run. 

3.2. Initial conditions 

Regarding the initial conditions, we concentrate on creating a set 
of identical clusters in all cosmologies studied. Thus, we started 
from a set of 17 clusters in the ACDM cosmology and adapted 
their initial conditions to the different dark-energy cosmologies 
investigated. 



3.2.1. The original set of initial conditions 

For the present analysis, we use cluster models obtained us- 
ing the so called re-simulation technique. It consists of re- 
simulating at higher resolution regions selected from an existing 
large-scale cosmological simulation. For this work, we selected 
as a parent simulat ion an A^-body run with 512^ particles in a 
boxof479/!-'kpc jYoshida et al.l200Uljenkrns et alJ2001l) . Its 
background cosmological model is spatially flat with Hm.o = 0.3 
and i^A = 0.7 at the final epoch, identified with redshift zero. 
The Hubble constant is /; = 0.7 in units of lOOkms^' Mpc^\ 
and the power spectrum normalisation was set to ag — 0.9. The 
particle mass was 6.8 x IO'^'/i^'Mq, and the gravitational soft- 
ening was chosen as 30/i^'kpc. From the output of this sim- 
ulation at z = 0, we randomly selected ten spherical regions of 
radius between 5 and 10/z^' Mpc, each containing either one or 
a pair of dark matter haloes, with mass larger than 3 x lO'^^M©. 
The total number of cluster-sized objects turns out to be 17. 

We constructed new initial conditions for each of these re- 
gions using th e ZIC softw are package {Zoomed Initial Condi- 
tions; Tormen et al.ll 19971) . The procedure for this construction 
is as follows. The initial positions of the particles in the simula- 
tion run define a Lagrangian region. The initial density field in 
this Lagrangian region is resampled by placing a larger number 
of particles than were originally present; on average, the number 
of high-resolution (HR), dark-matter (DM) particles is 10^ for 
our simulations. In this way, the spatial and mass resolution can 
be increased at will. The mass resolution of the re-simulation 
ranges from 2 x 10^ to 6 X \tfh-^Mp, per DM particle so as 
to have each cluster consists of approximately the same number 
of particles. The gravitational softening is given by a 5/i^'kpc 
cubic spline smoothing for all HR particles. 

The number of particles outside the HR Lagrangian region 
was reduced by interpolating them onto a spherical grid cen- 
tred on the geometrical centre of the HR region. An angular 
resolution of the grid between 3 and 5 degrees in both angular 
directions produces on the order of 100,000 macro particles of 
varying mass and gravitational softening. Extensive testing has 
shown that this number is enough to guarantee an overall good 
description of the original tidal field on large scales. 

The distribution of high resolution particles in the new initial 
conditions samples all the fluctuations of the matter power spec- 
trum realization of the original cosmological run, plus a new and 
independent realization of high frequency fluctuations from the 
original spectrum; in this way the power spectrum is extended 
up to the Nyquist frequency of the new HR particle distribution. 

The ZIC package has been widely tested and was used to pro- 
duce initial conditions for ma ny resimulations at medium to ex- 
tremely high resolution (e.g. Tormen et al. "1997^. ISpringel et alJ 
.2001b: Yoshida et aL.2002''.Stoehr et al.,2002) . 

3.2.2. Modifying the initial conditions for dark-energy models 

Since the goal of this work is studying the evolution of the same 
set of clusters in different cosmologies, we decided to adapt the 
initial conditions from the ACDM cosmology to all the dark en- 
ergy cosmologies we wish to investigate. This is done in two 
steps. In the first step, the initial redshift is adapted such that 
the rms density fluctuation amplitude today is the same in all 
models despite their different dynamics. The initial redshift Zini 
for any model is thus implicitly determined by the ratio of linear 
growth factors D+{z), 
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Fig. 2. — One individual cluster is shown at different redshifts (z = 2 to z = in steps of 0.5 from bottom to top) in different cosmologies (columns 
as labelled) all normalised to Og = 0.9 today. The panels illustrate that the clusters in different cosmological models arose from the same initial 
conditions and thus appear morphologically similar, but subtle differences are visible in detail. 



In a second step, the peculiar velocities of the particles in their 
Lagrangian initial region have to be modified such as to reflect 
the change in initial redshift. At early times T, they can be de- 
scribed by the iZel'dovichl i 1 9721) approximation. 
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Table 1 . — Parameters used for adapting the initial conditions 
of the ACDM model to the other cosmological models used in 
this paper Models marked with superscript ^ have their nns 
fluctuation amplitude Og adapted to the COBE normalisation, 
see Table|2| 



where q is the Lagrangian particle position, and 'P{q) is the ve- 
locity potential in Lagrangian space. Thus, the new velocities 
can be calculated starting from the displacement field 'Vg(S>{q) 
and scaling them to the new initial redshift and the new cosmo- 
logical model we wish to adapt them to. The velocities are thus 
given by 
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Table[2summarises the initial redshifts and the velocity scalings 
used for the different cosmologies. 
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3.3. Post-processing 



We then fit the NFW profile 



We applie d a halo finde r and c onstructed the merger tree as de- 
scribed in lTormen et'al] (120031) . 



3.3.1. Finding the halos 

The halo finder adopts the spherical overdensity criterion to de- 
fine collapsed structures in the simulations. This is done for each 
snapshot of each resimulation, estimating the local dark mat- 
ter density at the position of each particle, P,;dm^ by calculating 
the distance dj u) to the tenth closest neighbour, and assuming 



Pi;DM 



We then sort the particles by density and take as 



centre of the first halo the position of the particle embedded into 
the highest density. Around this centre, we grow spherical shells 
of matter, recording the total mean overdensity inside the sphere 
as it decreases with increasing radius. We stop the growth and 
cut the halo when the overdensity first drops below 200 times the 
mean (as opposed to critical) background density, and denote 
the radius so defined as r2oo- The particles selected in this way 
belong to the same halo and are used to compute its virial prop- 
erties (mass, radius, etc.). We tag all halo particles as engaged 
in the list of sorted densities, and selected the centre of the next 
halo at the position of the densest available (unengaged) parti- 
cle. We continue in this manner until all particles are screened. 
We include in our catalogue only such halos which have at least 
n = 10 dark-matter particles within their virial radius. All other 
particles are considered field particles. 

We found that our 17 clusters contain on average Nv ~ 
200,000 dark matter particles within their virial radii. The cor- 
responding virial masses range between My 
1.7 X \Q^^h-^MQ. 



3.1 X lO'^ to 



Po 



(10) 



f Navarro et al., 1997) to the binned profiles obtained. The con- 
centration is then defined as usual, c — r2oo/'"s- Note, however, 
that Navarro et al. ( 1997) defined r2oo as enclosing a mean den- 
sity of 200 times the critical rather than the mean background 
density as we do here. 
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Fig. 3. — Halo-mass histograms in the ACDM model for four of the 
eight equidistant redshift bins between redshifts zero and three intro- 
duced in Sect. |4] below. The distributions in the missing bins fall in 
between those in the adjacent bins and are left out for clarity. As halos 
merge and new halos form, the histograms shift towards higher masses. 



3.3.2. Constructing merger trees 

For each of these cluster-sized halos we build a merging history 
tree using the halo catalogs at all time outputs. Starting with 
a halo at any given z, we define its progenitors at the previous 
output z + Az to be all haloes containing at least one particle that 
by z will belong to the first halo. We call the main progenitor at 
z + Az the one giving the largest mass contribution to the halo at 
z. 



3.3.3. Fitting halo concentrations 

For determining halo concentrations and their change with time, 
we only took the main progenitor of each of the 17 massive ha- 
los into account. Using the centres from the halo finder, we 
construct radial profiles by binning the particles within r2oo in 
logarithmic radial bins with at least 132 particles each. For low- 
mass halos, this number is reduced to be 1/20 of the total num- 
ber of particles within r2oo to enshure enough points for the fit- 
ting procedure. Our choice in nu mber of particles is smaller 
than suggested bv'Powe r et alJ i200 3) but still reasonable, as we 
are interested in the local density at the innermost point, and 
not in the enclosed over-density. Also the global fit does not 
strongly depend on the innermost data-points. We checked, that 
the largest contribution to the scatter of the concentration pa- 
rameter within a redshift bin comes from the dynamical state of 
the halo. The changes in the concentration parameters infeiTed 
from the fits when changing then number of particles within the 
innermost bin by a factor of five is still one order of magnitude 
smaller than the intrinsic scatter due to the dynamical state of 
the halo. 



4. Halo concentrations 

Given the numerically determined halo concentrations, we can 
now proceed to compare them with expectations from analytic 
algorithms. Three such algorithms have been proposed: 

• iNavarro et al.l J 19971) assign a halo formation redshift Zcoii 
requiring that half of the final halo mass M be contained 
in progenitors with masses of at least /nfw M. The halo 
density scale po is then assumed to be Cnfw times the 
mean cosmic density at Zcoii- Fitting to their numerically 
determined halo concentrations, they recommend setting 
/nfw = 0.01 and Cnfw = 3 x 10^. 

• Bullock etal. (2001) noticed that the algorithm proposed 
by N avanxLet. al. ( 1997) . although well describing halo 
concentrations at redshift zero, under-predicts their de- 
crease with increasing redshift. They define the collapse 
redshift Zcoii by requiring that the non-linear mass scale at 
Zcoii be a fraction /b of the final halo mass, and assume that 
the halo concentration is a factor Cb times the ratio of scale 
factors at the redshift when the halo is identified and the 
collapse redshift. They recommend setting /b — 0.01 and 
Cb=4. 

• lEke et all ( 120011) define the collapse redshift using the 
power spectrum. They require that a suitably defined am- 
plitude of the dark-matter power spectrum at the mass scale 
M, linearly evolved to Zcoii, be equal to a constant Cens, and 
find Cens = 1/28. 



The algorithms by iNavarro et alJ ( H997l| ) and iBullock et alJ 
i200H ) involve two parameters, a mass fraction /nfw.b and a 
constant factor Cnfw.b, while Eke et al. (2001) introduce a sin- 
gle parameter only. 

Before we can apply these algorithms to our simulated ha- 
los, we have to introduce a common definition of virial radii and 
masses. While Navarro et al. define the virial radius r2oo 

as enclosing an overdensity of 200 ti mes the crit ical density of 
the Universe, Bullock et al. ( 2001 ) and iEkeetaP (5001) refer to 
the mean density instead. Since the mean is lower than the crit- 
ical density, r?on is l arger using the latter definition. Moreover, 
IBullock et al.l (l200Tt) use the mean overdensity within virialised 
spherical halos instead of the factor 200. 

In our numerical simulations, we use the radius enclosing 
200 times the mean background density as virial radius, which 
we call rjon th roughout. We adapt the numerical algorithms by 
INavarro et alj ll997). Bullock et al.. (2001 ) and Eke et al..(,2001) 
such that they use the same definition of raoo- This requires iter- 
ative solutions because converting masses from one definition of 
the virial radius to another requires the density profile, and thus 
concentrations, to be known. 

Given the masses and redshifts of all individual halos, we 
compute the analytically expected concentrations for each halo 
according to the three algorithms listed above. We thus have for 
each halo four concentration values, viz. the three analytic ex- 
pectations and the value obtained by fitting the NEW profile to 
the numerically determined density profile. We then bin halos by 
redshift into eight bins between redshifts three and zero. Each 
redshift bin thus contains halos with a range of masses. In each 
bin, we determine the median and the 33- and 68-percentiles of 
the distribution of the four types of concentration. The left panel 
of Eig.|4]shows results for the ACDM model as an example. 

The solid curve shows the numerically determined halo con- 
centrations, with the error bars giving the 33- and 68-percentiles 
centred on the median within each bin. The other curves show 
results obtained from the three algorithms as indicated in the 
figure. Error bars attached to those reflect the mass ranges of 
halos with in the respective redshift bins. While the algorithm by 
| Eke et all (.2001,.) describes the numerically fitted halo concen- 
trations very well within the error bars , concentration s obtained 
according to ' Navarro eTal] ( Il997h and lBuUock et alJ ( 1200 IB fall 
consistently above and below the numerical concentrations, re- 
spectively. The curve for the NEW algorithm reflects the earlier 
finding that it reproduces halo concentrations well at redshift 
zero, but under-predicts their decrease wit h increasing redshift. 
Concentrations according to .Bullock et alJ f2001) are somewhat 
too small in all redshift bins. 

The error bars on the simulated curve are larger than those 
on the analytically determined curves because of halo mergers, 
which change halo masses much more abruptly than halo con- 
centrations. 

Since all three algorithms have one or two free parameters, 
we now investigate whether their agreement with the numeri- 
cal results can be improved modifying the parameters. We thus 
define a measure for the quadratic deviation of analytical from 
numerical halo concentrations, 



X (^analytic ^numerical) 



halos 



(11) 



and minimise varying the parameter(s). Eigure |5l shows 
contours obtained f or the ACDM model from the a lgorithms by 
INavarro et a? (T997) (left panel) and Buflock et al. (2001) (right 
panel). Both algorithms have two parameters, a mass fraction 
and a factor, which are shown along the horizontal and vertical 
axes, respectively. 



Both contour plots share the feature of having a shallow, de- 
generate valley along which the mass fraction changes by or- 
ders of magnitude, while the factor changes very little. Approx- 
imately, these valleys follow the relations indicated by dashed 
lines in the two panels of Eig. [S] This near-insensitivity to the 
exact mass fraction allows fixing it at /nfw,b = 0.001, say, and 
then determining the minimum along the perpendicular axis 
only. 

Eigure |6l shows cuts through the contour plots in Eig.|5l and 
also accordin g to the one-parameter algorithm proposed by 
lEke et all ( 120011) . The abscissa shows the factors Cnfw.b di- 
vided by their values at the minima for fixed /nfw,b = 
0.001, i.e. Cnfw = 1500 and Cb ~ 3.5, while Cens is di- 
vided by the originally proposed value, i.e. Cens = 1 /28. The 
curves show pronounced minima for the algorithms proposed by 
Bullock et al. (2001) and Eke et al. (2001), and a very shallow 
minimum for the iNavarro et alJ(ll997i) algorithm. 
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Fig. 6. — Cuts through the % contours at constant mass fraction for 
three different halo-concentration prescriptions. 

Using these parameters, we now return to the comparing 
the numerically determined and analytically expected halo con- 
centrations and find the results shown in the right panel of 
Eig.m The agreern ent between the numerical results and the 
IBulloc k et al.' ('2001') concentrations has improved substantially, 
while the Navarro et al. ( 1997) algorithm still predicts too high 
halo concentrations at moderate and high redshifts. 

We repeated this procedure for all cosmological models used 
and found that the algorithm by Bullock et al. ( 2001 ) reproduces 
the numerical concentrations well in all of them if its parame- 
ters are m odified to fn = 0.001 and Cb = 3.5, while the algo- 
rithm by Eke et al. (2001) performs very well throughout with 
the factor Cens = 1/28 which was originally proposed. The 
agreement achieved for dark-energy models is illustrated for the 
Ratra-Peebles and SUGRA models in the left and right panels of 
Fig-El respectively. 

Our first conclusion is thus that the halo-concentration al- 
gorithms proposed by Bullock et al. (2001) performs very well 
also in dark-energy cosmologies, provided its two parameters 
are modified, while the algorithm by E ke et al. (.200 1) does not 
require any adaptati on. The algorithm originally proposed by 
INavarro et alJJl997l) has the same weakness in dark-energy as in 
ACDM models of under-predicting the redshift evolution. This 
can be remedied to some degree, but not removed, by modifying 
its two parameters. 

These algorithms have the advantage of relating halo proper- 
ties to the physical mechanism of halo formation. The gentle 



ACDM 



ACDM 




2 I , 1 , , 1 1 2 I ' ' ' ' ' 1 

0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

<z> <z> 

Fig. 4. — Concentrations of halos formed in the ACDM model binned by redsiiift before (left panel) and after (right panel) determining halo- 
prescription parameters. While the ENS algorithm fits the measured halo concentrations very well, concentrations according to the Navarro et 
al. and Bullock et al. algorithms deviate significantly. Changing the algorithm parameters yields very good agreement also for the Bullock et 
al. algorithm, while the Navarro et al. algorithm has too shallow redshift evolution. The right panel also shows the good quality of the simple 
polynomial fit. The error bars indicate the approximate l-o range of halo concentrations within the redshift bins. Concentrations are defined 
relative to the radius r2oo of spheres around halo centres containing 200 times the mean background density. 




Fig. 5. — contours in the space of parameters contained in the halo-concentration prescriptions by Navarro et al. (left panel) and Bullock et al. 
(right panel). Both contour plots show extended valleys in the surface with very flat bottoms. The dotted curves indicate approximate relations 
for the directions of these valleys. 



changes of halo concentrations with mass and redshift suggest, 
however, that they can be described in a much more simplified 
fashion by a fitting formula. Since concentrations of halos with 
fixed mass increase with decreasing redshift as the background 
density decreases, we attempt fitting the two-parameter func- 
tional relationship 

to the numerically determined halo concentrations. Best-fitting 
values for cq and a are listed for all eight cosmological mod- 
els in Tab. 121 The uncertainties of co and a given there are l-o 
bootstrap errors obtained from 1024 bootstrap halo samples. Al- 
though Eq. ( I12t does not reflect the physics of halo formation, it 
is useful for conveniently summarising the numerical results. 

The agreement between the numerically determined concen- 
trations and those obtained from the fitting formula il2\ are il- 
lustrated for ACDM in the right panel of Fig. @] and for two 



Table 2. — Parameters cq and a for the different cosmological 
models 



model 


08 


CO 


a 


ACDM 


0.90 


9.59 ±0.07 


-0.102 ±0.004 


Ratra-Peebles 


0.90 


10.20 ±0.07 


-0.094 ±0.005 


Ratra-Peebles 


0.82 


9.30 ±0.06 


-0.108 ±0.005 


SUGRA 


0.90 


11.15±0.09 


-0.094 ±0.006 


SUGRA 


0.76 


9.46 ±0.07 


-0.099 ±0.005 


OCDM 


0.90 


14.29±0.13 


-0.089 ±0.006 


w = -0.6 


0.90 


11. 32 ±0.09 


-0.092 ±0.005 


w = -0.6 


0.86 


10.44 ±0.08 


-0.066 ±0.005 



Ratra-Peebles 



SUGRA 
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Fig. 7. — Examples for the agreement between numerically simulated and analytically expected halo concentrations for two dark-energy models, 
Ratra-Peebles (left panel) and SUGRA (right panel). The agreement is very good, except for the Navarro et al. prescription whose redshift 
evolution is too shallow. 



dark-energy models in both panels of Fig. While the con- 
stant factor Co changes quite appreciably across cosmological 
models, the exponent a w —0.1 is approximately constant for 
all models, except perhaps for the dark-ener gy model wi t h con- 
stant w — —0.6 and reduced Og. Recently. IZhao et all J2003I) 
found that the dependence of halo concentrations on halo mass 
becomes shallower with increasing redshift. Since we are fo- 
cussing on how halo properties yary across cosmologies, our 
halo sample is cuiTently too small for confirming this result. 

Remarkably, halo concentrations are distributed about their 
mean yalues c in a way which is yirtually independent of the 
cosmological model. Figure [S] shows the distribution of c/c for 
all eight cosmologies, where separate curyes are shown for c 
determined according to the algorithms byBullock et al. ( 2001 ), 
lEke et alJ J2001I) and the fitting formula O. 
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Fig. 8. — Normalised halo-concentration distributions across all cos- 
mological models for three halo-concentration algorithms. The curves 
are very close to log-normal distributions whose standard deviations are 
almost independent of the cosmological model. 



These curyes closely resemble log-normal distributions, as 
found earlier in ACDM and OCDM cosmologies (e.g. Jina 
I200flt iBullock et al.ll200ll) . We thus attempt fitting them with 



distributions of the form 



p{c) = A exp 



(Inc — Inc — In/j) 



21 



2a2 



(13) 



The amplitude A is fixed by normalisation, and the parameters 
/J and o are determined from fitting to the measured distribu- 
tions, using c as obtained from the different algorithms and the 
polynomial of WIS . 

Any deyiation of /j from unity indicates that the respectiye dis- 
tribution is not centred on c. Since the algorithms and Eq. ( I12t 
were adapted such as to optimise the agreement between c and 
the measured mean yalues, /j is determined for cross-checking 
only. Table |3| shows the results for the three algorithms and the 
poly nomi al (Il2t . 

The table shows that the mean u scatters slightly about unity 
for all algorithms except for lNayarro et al.l ( il997.) . for which it 
is consistently below unity. This reflects the well-known inabil- 
ity of that algorithm to model the redshift eyolution correctly: 
It oyer-predicts concentrations c at redshifts aboye unity, and 
thus shifts the distribution of c/c systematically to yalues below 
unity. Perhaps more surprising is that the standard deyiation of 
the log-normal distribution is quite independent of the cosmo- 
logical model and of the algorithm used for describing the halo 
concentrations, a « 0.22. This result agrees well with the scatter 
found in a variety of CDM models with and without cosmologi- 
cal constant Llin g- .2000: .Bullock et al.-.200D . 

5. Dependence on w{z) 

In this section, we focus on the models with fixed Og — 0.9, 
in order to study the dependence of the halo concentration on 
the behaviour of vv(z). As we saw in the previous section, 
Co for fixed normalisation ag increases for the following order 
of cosmologies: ACDM, RP, SUGRA, DECDM, and OCDM. 
This is not accidental, but expresses one of our major results. 
As we stressed in the introduction, the different candidates for 
explaining the dark energy and cosmological acceleration are 
different in the redshift dependence of their equation of state. 
The models we study here, in particular the RP and SUGRA 
quintessence scenarios, well represent this issue, having the 
same wq but a markedly different w(z), as is obvious from Fig.[l] 
Thus, it is worthwhile studying the difference between these 
models in a phenomenological context. In Fig. |9l we plot the 



Table 3 . — Means and standard deviations of log-normal distributions fitted to normalised halo concentration distributions for four 
different halo-concentration prescriptions in the eight different cosmological models 



model 




mean 






a 








NFW 


B 


ENS 


poly 


NFW 


B 


ENS 


poly 


ACDM 


0.91 


1.04 


1.05 


1.00 


0.23 


0.22 


0.22 


0.21 


Ratra-Peebles 


0.93 


1.01 


1.03 


1.00 


0.23 


0.22 


0.22 


0.21 


Ratra-Peebles 


0.90 


1.00 


1.02 


1.00 


0.23 


0.22 


0.22 


0.22 


SUGRA 


0.95 


1.02 


0.98 


1.00 


0.24 


0.23 


0.23 


0.22 


SUGRA 


0.92 


1.02 


0.98 


1.00 


0.25 


0.24 


0.24 


0.23 


OCDM 


0.99 


1.10 


0.91 


1.00 


0.25 


0.24 


0.23 


0.23 


w = -0.6 


0.95 


1.00 


1.01 


1.00 


0.23 


0.23 


0.22 


0.22 


w = -0.6 


0.92 


0.99 


1.00 


1.00 


0.24 


0.23 


0.23 


0.22 


combined 


0.93 


1.02 


1.00 


1.00 


0.24 


0.23 


0.23 


0.22 



Table 4. — Concentration and perturbation asymptotic growth 
factor in the different cosmologies. 



cosmology 


CO 




D+(zcon) 


^ACDM 

cn 




L>+ (Zcoll ) 


ACDM 


1 


1 


1 


RP 


1.06±0.011 


1.07 


1.07 ±0.0003 


SUGRA 


1.16±0.013 


1.22 


1.19 ±0.004 


DECDM 


1.18±0.013 


1.19 


1.18±0.001 


OCDM 


1.49 ±0.018 


1.64 


1.61 ±0.01 



linear perturbation growth factor divided by the scale factor, 
g{z) = D^{z)/a{z), for the cosmologies studied here. is the 
growing solution of the general linear perturbation equation 



5 + 2-5-47i:Gp5 = 

a 



(14) 



which can, for w = — 1 and w = —1/3, be integrated once to 
yield °^ da^^ Jq d^^da . We normalise to unity at present. 

The asymptotic behaviour D+ a for a ^ means that the 
growth factor at early times in all cosmologies converges to the 
behaviour in a flat CDM Einstein-de Sitter case, with Ho = 1. 
On the other hand, the growth factor normalised by the scale 
factor, g{z), behaves differently for the different cosmologies. 
Most importantly, it turns out to reflect the behaviour of cq in 
the different cosmologies studied here above. 

The curves in Fig.|9]display g{z) for the different cosmologi- 
cal models, normalised to unity at present: the higher g{z) is at 
a given redshift z, the higher was the perturbation amplitude at 
that epoch, and the earlier was the epoch when structures formed 
(Bartelmann et al., 2002). Since at high redshifts all cosmolo- 
gies approach the flat CDM Einstein-de Sitter limit in which 
oc a, the structures acquire their fluctuation amplitude at red- 
shifts determined by the ratio of the asymptotic values of g{z) as 
shown in Fig.|9l 

It is thus natural to expect that the concentration parameters cq 
of our halos are affected similarly in the different cosmologies, 
since they parametrise the central density contrast of the halos. 
This is indeed what we find, as Tab.lSlshows. 

The table has four columns. The first abbreviates the cosmo- 
logical model. The second gives the cq parameters relative to 
^ACDM j.jjg /^(3j)]y[ model, with error bars obtained from the 
bootstrap errors of cq as given in Tab.|2l The third column shows 
the linear growth factors at infinity relative to the ACDM model, 
and the fourth column gives the same ratio taken at the average 



collapse redshifts of our numerically simulated halos. Error bars 
on the values in the fourth column are due to the mass range of 
halos which implies a range of collapse redshifts. 

The collapse red shifts were obtained using the prescription 
by Eke et al. ( 2001 ) because their algorithm for computing halo 
concentrations turned out to reflect our numerical results best. 
We draw three main conclusions from the results in Tab. |3 
First, the concentration parameters cq relative to ACDM are 
quite close to the linear growth factors at high redshift relative to 
ACDM, although the match is not perfect. We may be affected 
by cosmic variance because the number of halos per cosmology 
is relatively small. Second, using different definitions of the col- 
lapse redshift yields different ratios of the linear growth factor 
In particular, such definitions according to which the collapse 
redshifts are small fail in reproducing the trend in cq with dif- 
ferent cosmological models. Third, since the asymptotic ratio 
of growth factors towards infinite redshift is rather close to the 
ratios of the cq parameters, the final halo concentrations are ap- 
parently determined at very high redshift already even though 
the halos are quite massive and thus completely assembled late 
in cosmic history. 

We thus propose to interpret the cq parameter as composed 
of a factor valid for to the ACDM cosmology, multiplied by a 
correction which takes into account the asymptotic behaviour 
of the linear growth factors at high redshift for the given dark- 
energy cosmology. 



CO 



ACDM 
• "^0 



D 



D+{Z cou) 

(Zcoll) 



ACDM 



(15) 



This result is particularly important since it demonstrates that 
quintessence cosmologies having the same equation of state at 
present but a different redshift behaviour can produce relevant 
and predictable differences in the central regions of clusters. 

6. Summary 

We used cosmological numerical simulations to study the con- 
centration parameters of massive dark-matter halos in cosmo- 
logical models with dark energy, and compared them to results 
obtained in cosmological-constant (ACDM) and low-density, 
open (OCDM) models. 

We investigate three different dark-energy cosmologies. One 
has a constant ratio w — —0.6 between pressure and density of 
the dark energy, the other two have e quations of state which 
change over time. The Ratra-Peebles (Peebles & Ratra, 2002) 
model describes the dark-energy scalar field with an inverse 
power-law potential, the SUGRA (Brax & Martin, 2000) model 
modifies the power-law potential with an exponential factor 
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• Algorithms for predicting halo con centrations from thei r 
mass and redshift as proposed by 'Navar ro et alJ il997ft . 
[Bullock et al. (2001) and Eke et al. (2001) perform differ- 
ently in reproducing the concentrations measured in the 
simulations. The algorithms by Navarro et al. ( 1997) and 
jBullock et al. (2001) have two parameters each, a fraction 
of the final halo mass required in defining the collapse red- 
shift, and a factor relating core halo densities to the (mean 
or critical) background density at the collapse time. While 
the algorithm by Bullock et al. ( 2001 ) can be adapted mod- 
if ying their ori ginally proposed parameters, the algorithm 
by Navarro et aL (.1997) consistently predicts a too shallow 
redshift evolution of halo concentrations. The algorithm by 
lEke et alj (120011) succeeds very well in reproducing halo 
concentrations in all cosmologies studied without adapta- 
tion. 



Fig. 9. — Linear density perturbation growth dynamics in the different 
cosmologies considered. 



Dark energy with non-constant density, and more generally 
with non-constant equation of state, modifies terms in Fried- 
mann's equation. We modified and extended the proprietary 
GADGET-2 code accordingly, which is an extension of the pub- 
lically available GADGET code (Springel et al, 2001c). 

Cluster-sized halos were selected from a large cosmological 
simulation volume and re-simulated at substantially higher reso- 
lution. Since the large-scale simulation was originally prepared 
for a cosmological-constant (ACDM) cosmology, the initial con- 
ditions used for re-simulation needed to be carefully adapted in 
two ways to the other cosmological models. First, microwave- 
background (CMB) observations on large angular scales con- 
strain the amplitude of dark-matter fluctuations, conventionally 
expressed by the rms fluctuation level Og within spheres of 
8/i^'Mpc radius. Due to the enhanced dynamics of the grav- 
itational potential in dark-energy cosmologies, the integrated 
Sachs-Wolfe effect increases and thus the amount of primordial 
fluctuations on large scales needs to be reduced in order to re- 
main compatible with observed CMB fluctuations on large an- 
gular scales. This reduction of Og requires the redshift of the 
initial conditions to be adapted. Second, the particle velocities 
in the initial data need to be scaled according to the modified 
cosmological dynamics in dark-energy cosmologies. 

It should be noted that the correction of Og tries to adapt the 
density-fluctuation level in the simulations to CMB observations 
on large angular scales, while we are interested in the formation 
of cluster-sized halos on much smaller scales. How cluster-sized 
fluctuations are related to large-scale fluctuations depends sen- 
sitively on the primordial slope of the dark-matter power spec- 
trum, which we assume to unity here. 

We ran simulations for two sets of normalisations. As Og = 
0.9 for the original ACDM model, we first ran simulations with 
the same ag for the OCDM and the three dark-energy models 
in order to aiTive at cluster samples which could most directly 
be compared to each other We then ran three additional sets of 
simulations after adapting the normalisation of the dark-energy 
models to account for the enhanced integrated Sachs-Wolfe ef- 
fect. There are thus in total eight sets of simulations. Each of 
those contained 17 massive halos re-simulated at high resolu- 
tion. 

Concentrations for these halos were determined by fitting the 
NEW density profile at 78 redshifts between z = 5.8 and z = 0. 
Comparing these fitted values to analytical expectations on halo 
concentrations, we found the following results: 



• In all cosmologies studied, the halo concentrations permit 
the simple polynomial fit of Eq. Mil . While the mass de- 
pendence is shallow and quite independent of cosmology, 
the normalisation cq of the concentrations depends quite 
sensitively on the cosmological model. Halo concentra- 
tions are higher in models with time-varying dark energy 
than in cosmological-constant models. These findings con- 
firm numer ically what we had susp ected in an earlier ana- 
lytic study dBartelmann et alll2002h . 

• Once normalised to their analytically expected mean val- 
ues, halo concentrations closely follow a log-normal distri- 
bution with a standard deviation of a « 0.22, almost inde- 
pendent of the cosmological model. 

• The change of the concentration normalisation cq with cos- 
mology relative to the ACDM model is well represented by 
the respective ratio of linear growth factors at halo-collapse 
redshifts, e.g. as defined in the halo-concentration algo- 
rithm by Eke et al. (2001). Since this definition yields fairly 
high redshifts, this result seems to indicate that halo prop- 
erties are defined very early in cosmic history, in any case 
much earlier than the relatively late cosmological epoch 
when cluster-sized halos accrete most of their mass. 

Similar findings on the log-normal distribution of halo con- 
centrations and its standard deviation have be achieved before 
for CDM models with and without cosmological constant (see, 
e.g., Jing 2000; Bullock et al. 2001). Our study extends these 
investigations to dark-energy models. Our two most important 
conclusions are that (1) mean halo concentrations and their dis- 
tribution about those mean values are as in low-density open 
and cosmological-constant models except that (2) halos are sys- 
tematically more concentrated in dark-energy models in a way 
closely reflecting the linear density-fluctuation amplitude at high 
redshifts. 
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